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We discuss a property of our recently introduced generalized canonical ensemble (J. Stat. Phys. 
119 (2005) 1283). We show that this ensemble can be used to transform metastable or unstable 
(nonequilibrium) states of the standard canonical ensemble into stable (equilibrium) states within 
the generalized canonical ensemble. Equilibrium calculations within the generalized canonical en- 
semble can thus be used to obtain information about nonequilibrium states in the canonical ensemble. 
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The calculation of the equilibrium properties or states 
of systems having a nonconcave microcanonical entropy 
function (entropy as a function of the mean energy) is 
plagued by a fundamental problem. Because of the non- 
concavity of the entropy function, these systems possess 
energy-dependent equilibrium states in the microcanoni- 
cal ensemble that have no equivalent whatsoever within 
the set of temperature-dependent equilibrium states of 
the canonical ensemble 0,12, EJ 01 These so-called 
nonequivalent microcanonical states cannot, as a result, 
be calculated or predicted from the point of view of the 
canonical ensemble. To assess them, one must either re- 
sort to the microcanonical ensemble, with all the analyt- 
ical or numerical complications that this may imply, or 
else one must find a method with which one can study the 
nonequilibrium (i.e., metastable or unstable) states of the 
canonical ensemble. Indeed, it is known that nonequiva- 
lent microcanonical states correspond in general to meta- 
stable or unstable critical states of the canonical ensemble 
[H Q. El UDl < Hence, if one is able to calculate all the 
critical states of the canonical ensemble (equilibrium and 
nonequilibrium), then one has a complete knowledge of 
all the microcanonical equilibrium states, including the 
nonequivalent ones. 

Recently, we have followed this last line of thought to 
propose a generalization of the canonical ensemble which 
effectively transforms some of the metastable and unsta- 
ble critical states of the canonical ensemble into equilib- 
rium states of a generalized canonical ensemble [ll], [l2| . 
The generalized canonical ensemble is defined by the fol- 
lowing probability measure 



Z g {(3) 



,-aH{u)-ng{H{u)/n) 



where 



,(«) = J e-o^M-nrfffM/n)^ 
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represents a generalized canonical partition function. In 
these expressions, to denotes the microstate of an n-body 
system with Hamiltonian H(uS), a is an analog of an in- 
verse temperature, and g is some continuous function of 
the mean energy H(uS) / n which is as yet unspecified. The 
canonical ensemble is a special instance of this measure 
obtained, obviously, by choosing g = and a = (3. The 
case of quadratic functions g(u) = u 2 defines a statistical- 
mechanical ensemble known as the Gaussian ensemble 

mm. 

We have reported many properties of the generalized 
canonical ensemble in two recent papers 0, 0] — one 
containing rigorous results and their detailed proofs i 
and another one which is more "physical" in its presen- 
tation [T^ . The main point formulated in these papers 
is that the generalized canonical ensemble can be used, 
with suitable choices of the function g, to assess all the 
microcanonical equilibrium states of a system, including 
those not found at equilibrium in the canonical ensem- 
ble. In other words, for suitable choices of g, the gener- 
alized canonical ensemble can be made equivalent with 
the microcanonical ensemble in the thermodynamic limit. 
The conditions on g that ensure equivalence are stated 
in these papers both in terms of the microcanonical en- 
tropy function and in terms of a generalized canonical 
free energy function defined from Z g (a). 

Our aim in this short contribution is to discuss one in- 
teresting property of the generalized canonical ensemble 
that we have alluded to above, namely that the general- 
ized canonical ensemble can be used to transform meta- 
stable or unstable (nonequilibrium) states of the canon- 
ical ensemble into stable (equilibrium) states. We shall 
try to explain herein how this is possible and how this 
works in practice. 

The framework of our discussion is the same as in our 
previous papers 0, ^3 ■ We consider an n-body system 
described by some Hamiltonian function H(u>), with u) 
denoting the microstate of the system. In the canoni- 
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cal ensemble, that is, in a situation where the system is 
in contact with a fixed-temperature heat bath, all the 
equilibrium properties of the system are calculated using 
Gibbs's canonical measure 

p -PH(u) r 

w = ~zW m = J e ~ mu)du - (3) 

To calculate, for example, the equilibrium or thermody- 
namic-limit value of the mean Hamiltonian h(u>) — 
H(w)/n associated with a given value (3 of the inverse 
temperature, one first writes the probability measure for 
h, i.e., 

Pp(u) = I Pp(oj)duj = I 5(h(u) -u)Pp(uj)duj, 

J {uj:h(uj)—u} J 

(4) 

and then finds the global maximum of Pp(u) in the limit 
where n — > oo. The result of these steps is well-known: 
the global maximum of Pp(u), which corresponds again 
to the equilibrium value of h{w) in the canonical ensem- 
ble with inverse temperature [3, is given by the global 
minimum of the function 

Fp(u) = (3u- s{u) (5) 

evaluated at fixed [3. The function s{u) in this expres- 
sion is the microcanonical entropy function defined by 
the usual limit 

s(u) = lim — log / 8{h(uj) — u)dui. (6) 
n—>oo n J 

In the remainder, we shall denote the global minimum of 
Fp(u) at fixed (3 by up. 

This calculation of the equilibrium value up enables 
us to understand in a very simple manner how metasta- 
ble and unstable states can arise in the canonical ensem- 
ble. Also, we can understand with it why these states 
must arise in connection with nonequivalent ensembles 
and nonconcave entropies. The key points to observe 
here are the following (see also Fig.^J: 

(i) There can be mean-energy values u in the domain of 
Fp (it) that do not correspond to global minima of Fp (it) 
for any (3. In other words, the set of all possible values of 
Up, i.e., its range, can be a proper subset of the domain 
of Fp(u), which is the same as the domain of s(u). This 
happens when there is a first-order phase transition in 
the canonical ensemble; see Fig. ^ 

(ii) The domain of s(u) indicates which values of the 
mean energy are accessible to the microcanonical en- 
semble, while the range of up indicates which values of 
the mean energy are accessible to the canonical ensem- 
ble by varying (3. The previous point, therefore, implies 
the nonequivalcncc of the microcanonical and canonical 
ensemble — there are energy values accessible to the mi- 
crocanonical ensemble, but not to the canonical ensem- 
ble. 

(iii) Those values of u that are not global minima of 
Fp(u) for any (3 can be local minima or local maxima 



of Fp (it) for certain values of (3. The local minima cor- 
respond physically to metastable states of the canonical 
ensemble, while the local maxima correspond to unstable 
states. Note that the critical points of Fp(u) (minima, 
maxima or saddle-points) are given by Fp(u) = or, 
equivalently, s'(u) = (3 if s(it) is differentiable. 

(iv) In order to have local minima and local maxima 
of Fp(u), Fp(u) must be nonconvex, which implies that 
s(u) must be nonconcave. Accordingly, the nonconcavity 
of s(u) is a necessary condition for having metastable 
and unstable states in the canonical ensemble (it is also 
a sufficient one). 

These points have all been discussed at various lev- 
els of rigor in the literature; see, e.g., Ref. for a list 
of references on the subject of metastable and unstable 
states, and Ref. [l0| for more recent references on the 
connection with nonequivalent ensembles. For our pur- 
pose, what is important to keep in mind is that, for a 
given /3, Fp(u) may have local minima and local maxima 
in addition to one or more global minima correspond- 
ing to equilibrium states, and that all of these critical 
points of Fp(u) satisfy the equation s'(u) — (3 if s(u) is 
differentiable. Additionally, the stability of these critical 
points is determined by the sign of s"(it), provided that 
s(u) is twice differentiable. Indeed, by basic calculus, 
we have that if a point u satisfying F^(u) = is such 
that Fp(u) = —s"(u) > 0, then that point must be a 
minimum of Fp(u), either local of global. If the same 
point is such that Fp(u) — —s"(u) < 0, then it is a local 
maximum of Fp(u). 

Now, let us consider in more details the situation in 
which u is a microcanonically-allowed but canonically- 
unallowed value of h{ui) at equilibrium; that is, suppose 
that u is in the domain of s(u) but not in the range 
of up. As we have noted above, u can correspond to a 
metastable or unstable critical point of Fp(u), but never 
to a stable point of Fp(u), for the simple reason again 
that the stable points of Fp(u) correspond to the equilib- 
rium mean energies up in the canonical ensemble. How- 
ever, and this is the central point of this paper, u can 
correspond to an equilibrium mean energy of the general- 
ized canonical ensemble defined by Eg. P)). The analog of 
Fp(u) in the generalized canonical ensemble is the func- 
tion 

Fg fOI {u) = au + g(u)-s(u). (7) 

Therefore, the previous claim can be rephrased alterna- 
tively as following: there is a choice of function g and 
value a that make u a global minimum of F g ^ a (u). 

To demonstrate our claim, let us study the critical 
points of F g>a {u); these are defined by the condition 

F gta (u) = a + g'(u)-s'(u) = 0, (8) 

assuming that s(u) and g (it) are differentiable functions 
of u. Assuming further that these functions are twice dif- 
ferentiable, we can take the second derivative of F g ^ a {u) 
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FIG. 1: (a) Generic nonconcave micro-canonical entropy function, (b) The function Fp(u) constructed from s(u) shows, for 
some values of /3, a local minimum and a local maximum in addition to a global minimum. The global minimum defines the 
canonical equilibrium mean energy up, while the local min and local max correspond to metastable and unstable mean energies 
of the canonical ensemble, respectively, (c) Behavior of up as a function of (3. The nonconcavity of s(u) leads to a first-order 
phase transition in the canonical ensemble. This is reflected in Fp(u) by the fact that the global minima of this function never 
enter the range (ui,u h ). 



with respect to u to obtain 

F£ a (u)=g"(u)-s"(u). (9) 

From this equation, we readily see that the stability of the 
critical points of Fp(u) — F g= o. a= p(u) can be changed at 
will by choosing various functions g. Suppose, for exam- 
ple, that u is local maximum of Fp{u) satisfying s'(u) = (3 
and s"{u) > 0. Selecting g such that g"(u) > s"(u), 
and choosing a in the generalized canonical ensemble to 
be equal to the difference s'(u) — g'{ii), we have now 
that u must be a minimizer of F g ^ a {u), for then we have 
F' g a {u) = 0, and F' g ' a (u) > instead of Fg{u) < 0. This 
shows that a maximizer of Fp(u) can be changed to a 
minimizer of F gta (u). 

Whether or not u is a local or a global minimizer of 
F g , a {u) is undetermined at this point. To distinguish 
between the two types of minimizcrs, we need further 
notions of convex analysis that we shall not discuss here; 
see Ref. Q. However, it is possible to see with a specific 
example that the critical points of Fp{u) can be trans- 
formed into global minima of F g a (u) without any prob- 
lems. Indeed, consider a quadratic function g of the form 
g(u) — 7(u — ii) 2 , where u is a critical point of Fp(u). The 
function F g>a {u) of the generalized canonical ensemble 
has then the form 

F lt p(u) — au + "f(u — u) 2 — s(u). (10) 



Its first derivative evaluated at u is 

F^ {u) =a- s'(u) = 0, (11) 

so that the critical point u of Fp(u) is also a critical 
point of i ? 7 , Q =/3(w). Now to be sure that this critical 
point of F 7i /j(u) is a global and not just a local minimum 
of F 7i/ 3(m), we only have to choose a large, positive value 
for 7 that will make the term 7 (it — u) 2 dominate in 
-F 7)( g(u). Changing 7 in this case does not change the 
fact that u is a critical point of F lt<x [u), but it changes 
its "height" with respect to the other critical points of 
F 1<a (u). 

This transformation of the critical points of some func- 
tion to global minima of a different function involving an 
added quadratic term is well-known in the field of op- 
timization |lfj|. and has been described in the specific 
context of noncquivalcnt ensembles by Ellis, Haven and 
Turkington |l7j . Our work on the generalized canoni- 
cal ensemble can be seen as an extension of this previous 
work: it generalizes the quadratic-transformation trick to 
arbitrary penalty functions g in the way sketched above. 
At the statistical-mechanical level, the added function g 
leads, as we have seen, to the definition of a generalized 
canonical ensemble, and the basic property that g must 
possess in order to ensure that local minima or local max- 
ima of Fp(u) are properly transformed into global min- 
ima of F g ^ a {u) can be related physically to the existence 
of first-order transitions in the generalized canonical en- 
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semble [l2|. 

To understand this last point, recall that the metasta- 
ble and unstable mean energies of F g>ol {u) lie in a for- 
bidden range of mean-energy values for up, whose exis- 
tence defines physically a first-order phase transition in 
the canonical ensemble. By transforming these metasta- 
ble and unstable mean energies of the canonical ensemble 
to stable mean energies of a generalized ensemble, what 
one does, in effect, is to "shorten" or "inhibit" the range 
of forbidden mean energies in the generalized ensemble. 
That is, by choosing g, one seeks to make the jump of 
up seen in the canonical ensemble (Fig. (fl;) disappear in 
the generalized canonical ensemble. With this in mind, 
one can guess that the ultimate choice of g is one that 
makes the function s(u) — g(u) concave. In this case, 
Fg,a{u) must be concave for any value a, which implies 
that Fg i<x (u) can possess only one global minimizer for 
every a. Local minimizers or local maximizers are then 
not allowed, whereby it can be proved that one recovers 
full equivalence between the microcanonical and gener- 
alized canonical ensemble. The complete details as to 
which choices of g produce this effect can be found in 
Refs. see, e.g., Theorems 3 and 4 in Ref. 0], as 

well as Fig. 3 in that reference. 

To conclude this paper, we comment on four points: 

(i) Our discussion of metastability was confined to the 
level of the energy, but it can be generalized to describe 
metastability at the macrostate level as well. Many spin 
models are known, for example, to have magnetization 
states which are seen at equilibrium in the microcanoni- 
cal ensemble as a function of the energy, but not in the 
canonical ensemble as function of the temperature; see, 
e.g., Refs. H TTg| . Such nonequivalent magnetiza- 
tion states correspond to metastable or unstable magne- 
tization states in the canonical ensemble, and could, in 
theory, be transformed to equilibrium, stable states of a 
properly chosen generalized canonical ensemble. Work 
on this topic is ongoing ,20]. 

(ii) The generalized canonical ensemble (g ^ 0) can 
be interpreted physically as describing a system in ther- 



mal contact with a finite-size heat bath, as opposed to 
the canonical ensemble (g — 0) which describes systems 
in thermal contact with an infinite-size heat bath. For 
more details on this interpretation, the reader is referred 
to the original work of Challa and Hetherington which 
introduced the Gaussian ensemble 0, 0], as well as 
Refs. jHlll. 

(iii) The generalized canonical ensemble describes ex- 
tensive systems. In defining it, it is assumed that the 
Hamiltonian H scales like 0(n) in the thermodynamic 
limit, that In Z g (a) scales also like 0(n), and that the 
limit defining the microcanonical entropy in Eq.© ex- 
ists. 

(iv) Toral recently showed by direct calculation that 
the equilibrium properties of a mean-field spin model can 
be obtained by using a modified form of the partition 
function, which can be thought of to define a general- 
ized canonical ensemble [23|. His generalization of the 
partition function is related to our own generalization of 
the canonical ensemble (compare Eq. (4) in Ref. [23j with 
Eq. (PI of this paper) . Unfortunately, Toral does not pro- 
vide in his work any criteria for determining when the 
generalized canonical ensemble is actually useful, that is, 
when it is equivalent to the microcanonical ensemble or 
even the canonical ensemble. These criteria can be found 
in Refs. 00. 
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